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[i]  Radiation  belt  electrons  and  chorus  waves  are  an  outstanding  instance  of  the  important 
role  cyclotron  resonant  wave-particle  interactions  play  in  the  magnetosphere.  Chorus 
waves  are  particularly  complex,  often  occurring  with  large  amplitude,  narrowband  but 
drifting  frequency  and  fine  structure.  Nevertheless,  modeling  their  effect  on  radiation  belt 
electrons  with  bounce-averaged  broadband  quasi-linear  theory  seems  to  yield  reasonable 
results.  It  is  known  that  coherent  interactions  with  monochromatic  waves  can  cause 
particle  diffusion,  as  well  as  radically  different  phase  bunching  and  phase  trapping 
behavior.  Here  the  two  formulations  of  diffusion,  while  conceptually  different,  are  shown 
to  give  identical  diffusion  coefficients,  in  the  narrowband  limit  of  quasi-linear  theory.  It  is 
further  shown  that  suitably  averaging  the  monochromatic  diffusion  coefficients  over 
frequency  and  wave  normal  angle  parameters  reproduces  the  full  broadband  quasi-linear 
results.  This  may  account  for  the  rather  surprising  success  of  quasi-linear  theory  in 
modeling  radiation  belt  electrons  undergoing  diffusion  by  chorus  waves. 

Citation:  Albert,  J.  M.  (2010),  Diffusion  by  one  wave  and  by  many  waves,  J  Geophys.  Res  ,  115,  A00F05, 
doi.1 0.1 029/2 009J  A0 14732. 


1.  Introduction 

[2]  Cyclotron  resonant  wave-particlc  interactions  play  a 
key  role  in  both  the  acceleration  and  loss  of  radiation  belt 
electrons.  Chorus  waves,  in  particular,  are  believed  to  be 
key  to  both  the  energization  and  loss  of  energetic  electrons 
in  the  outer  zone  [ Chert  et  al.,  2007;  Horne ,  2007;  Bortnik 
and  Thorne ,  2007].  Chorus  waves  propagate  in  the  whis¬ 
tler  mode  and  arc  observed,  with  sufficient  time  resolution, 
to  be  coherent,  with  well-defined  frequencies  that  drift 
during  their  growth  to  large  amplitude  [Santolik  et  al.,  2003; 
Breneman  et  al. ,  2009].  The  wave  growth  is  intimately 
connected  to  the  linear  [Li  et  al ,  2008,  2009]  and  nonlinear 
[e.g.,  Nunn ,  1974,  Katoh  and  Omura ,  2007]  behavior  of 
resonant  electrons  with  energy  in  the  keV  range.  MeV  range 
electrons  are  also  subject  to  nonlinear  behavior  induced  by 
the  developed  waves,  but  their  motion  can  be  considered 
“parasitic,"  i.e.,  not  feeding  back  to  the  development  of  the 
waves. 

[3]  Coherent  cyclotron  resonant  interactions  of  test  elec¬ 
trons  with  individual  whistler  mode  waves  has  been  treated 
by  many  authors,  and  yields  three  distinctly  different  kinds 
of  particle  behavior,  namely  diffusion,  phase  bunching,  and 
phase  trapping.  Both  phase  bunching  (without  trapping) 
and  phase  trapping  are  favored  by  large  amplitude  waves 
and  low  inhomogeneity  of  the  background  magnetic  field;  a 
quantitative  criterion  has  been  developed  by  many  authors 
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[e.g.,  Inan  et  al ,  1978;  Albert ,  1993;  Omura  et  al .,  2008]. 
The  relevant  regime  also  depends  strongly  on  the  particle 
energy  and  pitch  angle,  so  all  three  types  of  behavior  may 
occur  under  the  same  conditions.  Albert  [1993,  2000, 
hereafter  Papers  I  and  II,  respectively]  derived  analytical 
expressions  for  the  changes  in  pitch  angle  and  energy  for  all 
three  types  of  motion,  using  a  Hamiltonian  formulation, 
though  frequency  drift  was  neglected.  Similar  considera¬ 
tions  also  apply  to  large  amplitude  electromagnetic  ion 
cyclotron  waves  [Albert  and  Bortnik ,  2009]  In  the  diffusive 
regime,  a  key  quantity  is  the  effective  interaction  time, 
which  is  controlled  by  how  long  (or  far)  the  particle  has  to 
move  in  the  varying  background  field  before  the  resonance 
condition  is  violated. 

[4]  The  large-scale  effects  of  chorus  waves  on  the  radia¬ 
tion  belts  have  also  been  modeled  using  quasi-linear  theory 
in  one,  two,  and  three  dimensions  (see  Albert  [2009]  for  a 
brief  review)  This  framework  assumes  a  eontinuum  of 
uncorrelated,  small  amplitude  waves,  with  wide  distribu¬ 
tions  in  frequency  and  wave  normal  angle,  in  a  eonstant 
background  magnetic  field.  Here,  the  diffusion  can  be 
considered  limited  by  the  relative  parallel  velocity  of  the 
particle  and  the  group  velocity  of  a  nearly  resonant  wave 
packet  [. Albert ,  2001].  The  resulting  local  pitch  angle  and 
energy  diffusion  coefficients  are  computed  locally  and  then 
bounce  averaged,  which  finally  introduces  variation  of  the 
background  magnetic  field  Recently,  the  expressions  for 
broadband  quasi-linear  diffusion  coefficients  were  expressed 
in  a  relatively  transparent  form  [ Albert ,  2005],  which  turned 
out  to  be  convenient  for  isolating  single  waves  within  the 
broad  frequency  and  wave  normal  angle  distributions.  Such 
single  waves,  suitably  chosen,  are  perhaps  surprisingly  well 
able  to  represent  the  entire  distributions,  leading  to  accurate 
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approximations  [Albert,  2007,  2008].  These  may  be  con¬ 
sidered  a  generalization  of  the  parallel  propagation  approx¬ 
imation  [Summers  et  ai ,  2007]. 

[5]  Thus  diffusion  emerges  from  both  quasi-linear  and 
nonlinear  treatments,  but  the  underlying  pictures  are  quite 
different.  Since  the  quasi-linear  diffusion  approach  seems  to 
model  the  actual  particle  behavior  fairly  well  [Albert,  2009], 
it  is  of  great  interest  to  relate  the  two  sets  of  diffusion 
coefficients  This  was  done  by  Albert  [2001],  working  with 
quasi-linear  expressions  for  whistler  mode  waves  in  the 
high-density,  low-frequency  limit  [Lyons  et  ai,  1972], 
which  invoked  considerable  simplifications  of  both  the 
whistler  mode  dispersion  relation  and  the  resonance  condi¬ 
tion.  It  was  concluded  that  the  narrowband  limit  of  the 
quasi-linear  pitch  angle  diffusion  coefficient  was  approxi¬ 
mately  equal  to  the  Hamiltonian-derived  pitch  angle  diffu¬ 
sion  coefficient  for  monochromatic  waves.  Here,  the 
comparison  of  the  two  analytical  frameworks  is  recon¬ 
sidered  in  much  greater  generality,  using  the  full  description 
of  stationary  cold  plasma  waves.  It  is  shown  that  the  nar¬ 
rowband  limit  of  bounce-averaged  quasi-linear  theory  and 
the  diffusive  regime  of  the  Hamiltonian  analysis  yield 
exactly  the  same  pitch  angle,  energy,  and  cross  diffusion 
coefficients  Furthermore,  averaging  the  monochromatic 
results  over  distributions  of  frequency  and  wave  normal 
angle,  which  statistically  models  a  sequence  of  resonant 
interactions  with  individual  waves,  recovers  exactly  the  full 
broadband  quasi-linear  diffusion  coefficients.  This  seems  a 
meaningful  step  toward  reconciling  the  behavior  expected 
from  coherent  dynamics,  in  the  diffusion  regime,  with  the 
apparent  utility  of  bounce-averaged  quasi-linear  theory  for 
modeling  of  radiation  belt  electrons. 

[6]  Any  possible  coupling  between  changes  in  cto  and  p 
with  changes  in  L  will  be  ignored.  This  is  usually  justified 
by  the  wide  separation  of  time  scales  associated  with  the 
first  two  adiabatic  invariants  compared  to  that  of  the  third, 
i.e ,  the  drift  period  compared  to  the  cyclotron  and  bounce 
periods.  Such  coupling,  which  would  lead  to  cross  diffusion 
terms  involving  D n(tL  and  DpL ,  has  only  been  considered 
occasionally,  usually  in  the  context  of  so-called  drift  shell 
splitting  [Roederer,  1970;  Schulz  and  Lanzerotti ,  1974], 
although  Brizard  and  Chan  [2004]  recently  formulated  the 
“full”  matrix  of  diffusion  coefficients  generated  by  an  ar¬ 
bitrary  wave  spectrum  in  axisymmetric  geometry.  The 
resulting  diffusion  equation  could  be  solved  numerically  by 
an  algorithm  based  on  stochastic  differential  equations  [Tao 
et  ai,  2008]  or  the  layer  method  described  by  Tao  et  ai 
[2009]. 

[7]  Section  2  exhibits  the  local  quasi-linear  diffusion 
coefficients  and  their  monochromatic  limit,  following  Albert 
[2007],  and  carries  out  the  bounce  average  following  Albert 
[2001],  leading  to  closed  form  expressions  with  no  re¬ 
maining  integrals.  Section  3  presents  the  diffusion  coeffi¬ 
cients  of  Albert  [1993,  2000]  resulting  from  coherent 
interactions  with  a  single,  monochromatic  wave,  which  are 
found  to  be  identical  to  the  final  results  of  section  2.  Section 
4  then  considers  the  coherent  diffusion  coefficients  suitably 
averaged  over  wave  frequency  and  wave  normal  angle 
parameters,  reproducing  the  fall  quasi-linear  expressions. 
Section  5  presents  some  numerical  examples  of  diffusion 
coefficients  for  a  model  of  nightside  chorus  waves,  calcu¬ 


lated  from  each  approach,  and  explicitly  demonstrates  their 
equivalence.  This  is  followed  by  a  brief  discussion. 


2.  Quasi-Linear  Diffusion  Coefficients 

[s]  The  condition  for  gyroresonance  between  a  particle 
and  a  wave  is 

u-kfV jjsft,,  n„=snilc/y,  (I) 

where  n  is  an  integer,  s  =  ±1  is  the  sign  of  the  charge  of  the 
particle,  Qc  =  | q\  Bimc  is  its  local  nonrclativistic  gyroffc- 
quency,  and  7  is  its  relativistic  factor.  The  local  pitch  angle 
of  the  particle  is  a,  the  index  of  refraction  is  p  =  kduj,  and 
the  wave  normal  angle  is  9.  The  underlying  mechanism  of 
quasi-linear  diffusion  can  be  thought  of  as  involving  con¬ 
tinuous  resonance:  even  as  the  particle  diffuses  in  a  and  7,  it 
is  always  able  to  find  an  instantaneously  resonant  wave 
within  the  u>  and  9  distributions. 

[9]  Albert  [2001]  considered  whistler  waves,  using  ex¬ 
pressions  based  on  the  approximations  uj/Qc  1  <£ 

[Lyons  et  ai,  1972;  Lyons,  1 974b],  but  here  any  cold  plasma 
mode  is  considered,  without  any  such  approximations. 

2.1.  Local  Expressions 

[10]  The  local  diffusion  coefficients  in  a  spectrum  of 
waves  were  given  by  Lyons  [1974a,  1974b],  as  derived  from 
the  Vlasov  equation  [Kennel  and  Engelmann ,  1966;  Lerche, 
1968],  although  it  can  also  be  obtained  by  considering 
motion  of  a  single  particle  acted  on  by  single  wave,  for  an 
interaction  time  related  to  the  wave  packet  bandwidth  [e.g., 
Albert,  2001].  In  either  case,  the  spatial  variation  of  the 
background  magnetic  field  and  all  other  parameters  is  ig¬ 
nored  for  the  local  calculation,  and  accounted  for  later  by 
bounce  averaging. 

[11]  The  derivation  is  fairly  involved  (see  also  the  pre¬ 
sentations  by  Walker  [1993]  and  Swanson  [1989]),  but  the 
results  for  pitch  angle  cx  and  momentum  p  can  be  expressed 
as 

^  7T  f  d^k  .  „  s 

2 J  ~ *||V|1  ~ n") 

;  I  Bit  I2  j$,,|2  (- sin2  a  +  n„/a;)2 

^  V  p2  P2  COS2  Q  ’ 

ff1,?  _  p  sin  a  cos  a  _  / D^\ 1 

-  sin2  a  +  n„/u/  /?„,  W,„/  P) 

D,m  has  dimensions  of  1//,  because  of  the  explicit  division 
by  p2  Bu  is  the  Fourier  transform  of  the  wave  magnetic  field 
taken  over  the  plasma  volume  V  (which  is  effectively  infi¬ 
nite),  and  J4>„|2,  as  given  by  equation  (9)  of  Lyons  [1974b], 
is  the  result  of  resonance  averaging  the  geometric  details  of 
the  particle  motion  in  the  electromagnetic  field  of  an  oblique 
plane  wave  The  ratios  of  the  diffusion  coefficients  were 
interpreted  by  Kennel  and  Engelmann  [1966]  in  terms  of 
single-wave  characteristics  of  a  quasi-linear  diffusion  op¬ 
erator,  and  were  farther  discussed  by  Lyons  [1974a]  and 
Summers  et  ai  [1998]. 

[12]  The  expressions  get  more  involved  after  transforming 
the  integration  variables  from  (Arj_,  k  )  to  (lj,  9 ),  and  mod- 
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eling  |Bk2/F  as  a  function  of  (lj,  6 ),  which  brings  in  nor¬ 
malization  integrals.  As  expressed  by  Albert  [2005],  and 
similarly  by  Glauert  and  Home  [2005],  the  resulting  form 
of  the  diffusion  coefficients  can  be  written  as  the  sum  over  n 
of  terms  Dn  given  by 


6  and  z.  As  the  0  distribution  is  narrowed,  ^rcs  becomes  a 
well-defined  function  of  z.  And  as  the  u  distribution  is 
narrowed,  B2^)  approaches  a  6  function  of  v.  Assuming  9m 
and  wm  are  compatible  with  resonance  at  some  location  zm, 
the  bounce  average  and  Gx  combine  to  give 


with 


Dn 


Q c  ^wavc 

72  B 2 


L 


sin  9d6  A„  G\  , 


(3) 


J—  Gi(wrcs(z,0m))  =*•  J 

-J 


—  f2c<5(u>TO(r)  -  uim) 
V||  ‘  \du)fdz\ 


(7) 


A„ 

G, 


C2 


r 


7 r  sec 6  (- sin2  a  + 

;|1  -(3w/a*,Vv:lr 

QcB2(lu) 


&Xg) 

fodff  sin  ff  g*,(&)  T 


=  M2 


M  4-  w 


(4) 


The  full  wave  intensity,  £wavc,  is  now  considered  to  be 
concentrated  at  the  single  pair  ( 6m ,  u;m). 

[16]  The  derivative  of  u)  is  evaluated  using  the  resonance 
condition,  and  it  is  important  to  note  that  k  is  a  function  of 
both  z  and  as  specified  by  the  dispersion  relation. 
Therefore  implicit  differentiation  of  the  resonance  condition 
gives 


(8) 


[13]  The  refractive  index  p  is  a  known  function  of  (*<;,  9) 
for  the  given  cold  plasma  wave  mode  [e.g.,  Stix ,  1962]. 
B2(uj)  describes  the  frequency  distribution  of  wave  power, 
and  is  nonzero  only  between  lower  and  upper  cutoffs,  u)LC  < 
uj  5  ^UC‘  Similarly,  the  distribution  of  wave  power  with 
wave  normal  angle  6  is  described  by  gjfi\  which  is  nonzero 
only  for  i9min  <  6  <  9miX.  Both  ^(J)  and  g^(t anO)  are  usually 
modeled  as  truncated  Gaussians,  peaked  at  um  and  9m , 
respectively.  The  quantities  G\  and  G2  are  explicitly  nor¬ 
malized  versions  of  B2(J)  and  gjfi),  and  are  discussed 
further  in  Appendix  A. 

2.2.  Narrowband  Limit 

[14]  As  shown  by  Albert  [2007,  2008],  the  integral  in 
equation  (3)  may  be  approximated  as  a  weighted  average, 
which  becomes  exact  as  gJJ))  becomes  narrowly  peaked.  In 
that  limit, 


Dn 


f L  ^wavc  ^ n 


(5) 


evaluated  at  some  resonant  pair  (u;,  6)  within  the  specified 
distributions.  For  the  purposes  of  Albert  [2007,  2008],  wLC 
and  uuc  were  used  to  find  6  ranges  containing  resonances, 
and  ZXU  was  approximated  using  representative  values  from 
within  these  ranges.  In  section  2.3,  equation  (5)  is  evaluated 
at  0mt  with  u  taken  to  be  the  corresponding  resonant  value  at 
each  location. 


[17]  The  factors  of  (An/T)/  dujldz\  containing  partial  deri¬ 
vatives  combine  and  simplify: 


1  du 

1 

i+-i 

p  di 

f  1 

dk 

V 

vll 

* 

11  dw 

uJ 

Putting  everything  together  gives 


2 T/,  B-  | V||  |v272  fi2  B  cos’  <>o 


:(-sin2"+v) 


9  U 


4-  fL) 


(10) 


where  Beq  and  ao  are  equatorial  values  but  all  other  quantities 
are  evaluated  at  the  resonance  location.  This  is  the  mono¬ 
chromatic  limit  of  the  bounce-averaged,  broadband  quasi- 
linear  diffusion  coefficient  for  each  n. 

[is]  The  bounce-averaged  coefficients  D^p  and  Ifpp  are 
derived  similarly,  and  in  the  monochromatic  limit  are  related 
to  D”^  by 


d: 


rr 

<V>riO 


p  sin  ao  cos  ao  B 
-  sin2  ct+VLn/uj  Beq 


D"n 

pp  _ 

D" 


(in 


for  each  n.  Albert  [2004]  discussed  the  role  of  these  ratios  in 
enforcing  the  condition  Dn<fXt  Dpp  >  D2<oP. 


2.3.  Bounce  Averaging 

[15]  The  bounce-averaged  diffusion  coefficient  for  the 
equatorial  pitch  angle,  ao,  is  given  by  the  sum  over  n  of 


[T 


^\daj 


(6) 


where  z  is  distance  along  the  magnetic  field  line  (and  is 
easily  converted  to  latitude).  In  equation  (3),  B2(u>)  is 
evaluated  at  the  resonant  frequency,  which  depends  on  both 


3.  Coherent  Interactions 

[19]  A  quite  different  scenario  is  that  of  a  particle  inter¬ 
acting  with  a  single  wave  in  a  spatially  varying  magnetic 
field,  so  that  the  resonance  condition  of  equation  (1)  is  only 
satisfied  at  discrete,  isolated  locations  through  which  the 
particle  passes.  As  mentioned,  analytical  estimates  of  the 
resulting  particle  motion  were  obtained  in  Papers  I  and  II 
For  large  amplitude  waves  and  small  background  inhomo- 
geneity,  nonlinear  behavior  (phase  bunching  and  phase 
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trapping)  can  occur,  but  here  the  opposite  limit  is  consid¬ 
ered,  which  leads  to  random  walks,  or  diffusion 

[20]  Papers  I  and  II  write  out  the  full  equations  of  motion 

in  Hamiltonian  form,  transform  to  gyroresonance  variables, 
expand  to  first  order  in  5wave/5,  and  appropriately  average 
away  nonresonant  terms.  For  n  £  0,  this  leads  to  two  con¬ 
stants  of  motion  which  can  be  used  to  reduce  the  number  of 
variables  to  a  single  action-angle  pair,  (/,  Q.  To  lowest  order, 
/is  proportional  to  the  familiar  first  adiadatic  invariant,  and 
f  is  the  usual  wave-particle  phase  which  is  stationary  at 
resonance.  The  evolution  equations  for  /  and  £  can  be 
expressed  in  terms  of  a  reduced  Hamiltonian,  K  =  K0(I ,  z)  + 
K\{Iy  z)  sinf,  with  z  playing  the  role  of  time  The  adiabatic 
motion  is  described  by  Kq>  while  K\  captures  the  effects  of 
the  resonant  wave.  For  n  —  0,  a  similar  treatment  yields  a 
reduced  Hamiltonian  M  —  A/0(Y,  z)  +  A/i(Y,  z)  sin£,  where 
T  =  72.  The  reduced  Hamiltonians  can  be  used  to  derive 
analytic  approximations  to  the  resonant  changes  in  the  adi¬ 
abatic  invariants  /  or  T.  An  “inhomogeneity  parameter”  R, 
proportional  to  (^5/dz)/Z?wave,  delineates  diffusion  from  the 
nonlinear  regimes  involving  phase  bunching  and/or  phase 
trapping.  Here  we  only  consider  the  case  \R\  1,  which 

indicates  diffusion. 

3.1.  Cyclotron  Resonance 

[21]  At  an  isolated  resonance  n  £  0,  according  to  Papers  I 
and  II, 


^02/&a/cos2(^  +  g/3-  (12) 

Again,  z  is  distance  along  the  field  line,  and  07  is  the  sign  of 
P2  KJdzdl  at  resonance.  Averaging  over  £rcs,  which  depends 
on  the  gyrophase  and  is  randomized  between  bounces, 
yields  1/2.  Papers  I  and  II  also  give  the  perturbation  Ham¬ 
iltonian  K\  in  terms  of  an ,  which  describes  the  wave  field 
components.  The  relation  between  Ku  0*,  and  noted  by 
Albert  [2001]  holds  for  general  cold  plasma  waves 


(16)  into  equation  (15)  yields  the  first  major  result  of  this 
paper:  the  coherent  interaction  versions  of  Z)”o nQi  D^o/n  and 
ITpp  work  out  to  be  exactly  the  same  as  in  equations  (10)  and 
(11)  for  the  narrowband  limit  of  the  bounce-averaged  quasi- 
linear  expressions. 

3.2.  Landau  Resonance 

[23]  For  the  special  case  n  =  0,  Paper  II  gives 

(6T)  - M'W^kd T  cos2(^  +  ^3'  <17) 

and  ct y  is  the  sign  of  T  at  resonance  Here  a0  is 

just  an  with  n  =  0,  but  now 

Ml=ain2  cos:0  =  4cos2^^^4>5.  (18) 


The  Hamiltonian  equation  of  motion  for  £  yields 


&M0  =  d  /rfA 
d:0T  d:\dz) 


c2  d  .  v 


(19) 


where  again  can  be  omitted  in  the  z  derivative. 
[24]  The  diffusion  coefficients  arc  now 


nr  0  _  «AT)2) 

2  Th 

m 1 

and  so  on.  Using 

da0 

2  ■> 
nrtr 

dp 

~df 

~  -  tan  «0  — y  1 

2 P2 

dr  ~~ 

(20) 


(21) 


from  Paper  II,  the  resulting  coherent  interaction  expressions 
for  and  Dp~°  again  agree  exactly  with  equations 

(10)  and  (11)  from  the  narrowband  limit  of  bounce-averaged 
quasi -linear  theory. 


K 


2 

1 


a*  »2  vj 

4  (Pu/mc)2  V2  v2  B- 


(13) 


The  Hamiltonian  equation  of  motion  for  £  yields 


cfKp  d_  f  d£\ 

dzdl  ~  dz\dz) 


c2  d  . 


■n. 


(14) 


where  u>  is  the  constant  frequency  of  the  single  wave  and 
can  be  omitted  in  the  z  derivative. 

[22]  Diffusion  coefficients  are  constructed  from 


4.  Average  Over  Wave  Distributions 

[25]  It  has  just  been  shown  that  the  monochromatic  limit 
of  bounce-averaged,  broadband  quasi -linear  theory  is  well 
behaved,  and  reduces  to  the  results  of  a  Hamiltonian  anal¬ 
ysis  of  a  resonant  interaction  with  a  single  wave  (in  the 
diffusive  regime).  Conversely,  pitch  angle  diffusion  by  a 
single,  coherent  wave  can  be  expressed  in  terms  of  the 
quantities  defined  for  quasi-linear  diffusion 


[F 

«0O(i 


^wavc  ^ O 

72  b 2  r 


Um)  (22) 


2 


dctQ  dp 

Wdi' 


(15) 


where  (  )  denotes  the  average  over  frcs.  From  Paper  II, 


for  either  n  ^  0  or  n  =  0.  Wc  now  consider  the  result  of  many 
coherent  interactions  with  individual  waves  all  with  ampli¬ 
tude  5wavc  but  with  frequency  and  wavenormal  angle  sta¬ 
tistically  distributed  according  to  B2{J)  and  gJQ). 

[26]  The  appropriate  average  is 


d(*o  _  —  sin2  (\+Qn/uj  Beq  m2cr  7  dp  m2c 2  7 

dl  sin *>o cos r*o  &  p2  sn'  dl  p  sn' 


'  J  |Bk  tfk 


(23) 


The  corresponding  ratio  dp/da q  is  closely  related  to  the  ra-  where  Dk  refers  to  the  single-wave  equation  (22).  The 
tios  in  equation  (11).  Substituting  equations  ( 1 2>— ( 14)  and  denominator  of  (23)  is  just  Z?£,ave.  Converting  from  to 
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Figure  1.  Local  quasi-linear  pitch  angle  diffusion  coeffi¬ 
cient  for  1  MeV  electrons  interacting  with  a  broadband  spec¬ 
trum  of  chorus  waves  at  L  ~  4.5  at  different  latitudes.  At  the 
equator,  the  spectrum  is  peaked  at  9m  =  0,  um/£le  =  0.35. 
Only  the  lowest  harmonic  (/?  =  —!)  term  is  shown. 


dudQ  in  the  numerator,  using  the  results  of  Appendix  A, 
gives 


(PLJ 


=  f 


dujd9D\ 


B2(u) 


sin  9gJ9)  T 


J  kfBHu)d'J  fdff  sin  O'  T 


(24) 


Then,  schematically, 

J  djj  J  d9  J  dz  6(um(z)  -  wm) 

=>  J  dz  J  d9  J  duj  6{ u>  —  ujm)=>  J  dz  j  d9 

which  yields 

- kit  (&)’ t %*  / 


This  is  the  second  major  result  of  this  paper:  the  coherent 
interaction  diffusion  coefficient,  suitably  averaged,  is  iden¬ 
tical  to  the  full  quasi-linear  result  given  by  equation  (6).  The 
analogous  relations  hold  for  (D?XoP)  and  (Dpp). 


5.  Numerical  Example 

[27]  For  illustration,  we  consider  the  model  of  Li  et  at. 
[2007]  for  nightside  chorus  during  a  magnetic  storm  main 
phase,  at  L  =  4.5  with  =  3.8  at  the  equator.  They 

computed  quasi-linear  diffusion  coefficients  for  waves  with 
#wavc  =  50  pT,  with  the  equatorial  frequency  distribution 
specified  by  =  0.35  ftc,  6uj  =  0.15  f^,  uLC  =  0.050*.,  and 
u> uc  =  0-65  Ot,.  The  waves  are  considered  present  only  for 
latitude  A  <  15°.  In  that  work  the  waves  were  all  taken  to 
propagate  with  0=0,  but  here,  following  Horne  et  al 


[2005]  and  Albert  [2008],  the  wave  normal  angle  distribu¬ 
tion  is  modeled  with  6m  =  0,69  =  30°,  0min  =  0,  and  0max  =  45°. 

[28]  Figure  1  shows  the  local  quasi-linear  pitch  angle 
diffusion  coefficients  for  1  MeV  electrons  for  several  values 
of  latitude,  calculated  from  equation  (3).  Only  contributions 
by  n  =  -1  are  shown.  For  each  wave  normal  angle  in  the 
distribution,  the  resonant  frequency  is  found;  if  both  lie 
within  the  model  distributions,  a  contribution  is  made  to  the 
diffusion  coefficient  integrals.  The  ‘usual*  quasi-linear 
results  [e.g  ,  Horne  et  al.,  2005;  Albert,  2005]  consist  of  just 
such  calculations,  converted  from  a  to  ao  and  bounce 
averaged,  as  in  equation  (6),  and  summed  over  n 

[29]  Figure  2  shows  equatorial  pitch  angle  diffusion 
coefficients  for  individual  waves  with  0  =  9m  =  0  and  var¬ 
ious  frequencies  between  ujlc  and  u>uc,  calculated  according 
to  the  coherent  formulation,  equations  (12)  and  (15),  with 
n  =  —  1.  Related  calculations  were  previously  presented  by 
Albert  [1993,  2000,  2002]  and  Albert  and  Bortnik  [2009], 
As  mentioned  in  section  3,  integration  along  the  field  line  is 
inherent  in  the  formulation.  The  curve  for  6  ~  9m,  to  —  ujm,  is 
emphasized  by  the  dashed  curve  Figure  3  is  similar,  but 
shows  the  results  holding  uj  =  ujm  fixed  and  varying  9  from 

$min  tO  0max- 

[30]  Figure  4  shows,  as  solid  curves,  the  quasi-linear 
diffusion  coefficients  after  carrying  out  the  bounce  averages 
of  the  local  results  illustrated  in  Figure  1.  The  sum  of 
contributions  from  n  =  - 1  and  n  =  +1  are  shown  in  the  top 
row,  and  just  n  =  0  is  shown  in  the  bottom  row.  Also  shown, 
as  red  squares,  are  the  results  of  numerically  averaging  the 
diffusion  coefficients  for  monochromatic  waves,  from 
Figures  2  and  3,  weighted  according  to  equation  (24).  It  is 
apparent  that,  allowing  for  numerical  accuracy,  the  com- 


D„oao  (day-1) 


ao 

Figure  2.  Equatorial  pitch  angle  diffusion  coefficient  for 
1  MeV  electrons  interacting  with  monochromatic  chorus 
waves  at  L  =  4.5,  treated  as  a  coherent  interaction.  Results 
are  shown  for  a  fixed  value  of  wave  normal  angle  and  sev¬ 
eral  fixed  values  of  frequency;  the  dashed  line  indicates  9  = 
6m  =  0,  uj  =  ujm  =  0.35f^  (at  the  equator).  Only  the  lowest 
harmonic  (n  =  -1)  term  is  shown. 
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D„oao  (day-1) 


Figure  3.  Same  as  Figure  2  but  showing  results  for  fixed 
frequency  and  several  values  of  wave  normal  angle 

putational  evaluations  verify  the  analytical  result  that  the 
two  formulations  are  the  same. 

6.  Summary  and  Discussion 

[31]  This  paper  has  investigated  the  relationship  between 
two  seemingly  different  formulations  of  wave-particlc 


interactions.  Generalizing  a  previous  study,  it  has  been 
shown  analytically  that  taking  the  narrowband  limit  of 
bounce-averaged,  broadband  quasi-linear  diffusion  coeffi¬ 
cients  agrees  exactly  with  the  diffusive  limit  of  coherent 
interactions  with  a  monochromatic  wave.  Moreover,  con¬ 
sidering  the  individual  waves  to  be  drawn  from  specified 
frequency  and  wavenormal  angle  distributions,  and  aver¬ 
aging  diffusion  coefficients  accordingly,  reproduces  the  full 
quasi-linear  expressions. 

[32]  It  has  been  a  puzzle  why  global  simulations  using 
quasi-linear  theory  [Li  et  ai,  2007;  Albert ,  2009]  are  at  least 
moderately  successful  in  reproducing  the  observed  effects  of 
chorus  waves,  which  upon  close  examination  are  discrete 
and  coherent  [Santolik  et  ai ,  2003].  Parameters  used  to 
model  chorus  waves  as  a  population  which  are  based  on 
wave  measurements  with  coarse  time  resolution  [Meredith 
et  ai ,  2003]  should  reflect  the  distribution  of  the  underly¬ 
ing  individual  waves.  As  just  shown,  multiple  interactions 
with  this  distribution  of  waves  will  be  well  described  sta¬ 
tistically  by  the  quasi-linear  approach,  as  long  as  the  indi¬ 
vidual  waves  are  not  large  enough  to  induce  nonlinear 
particle  behavior  [Cattell  et  ai ,  2008;  Cully  et  ai ,  2008]. 

[33]  It  should  be  noted  that  in  all  cases,  the  wave  para¬ 
meters  (amplitude,  frequency,  wave  normal  angle)  have  been 
treated  as  constant  during  each  individual  wave-particle 
interaction.  Although  the  quantities  can  vary  significantly, 
indeed,  frequency  drift  is  a  characteristic  feature  of  chorus 
waves,  the  duration  of  an  isolated  interaction  is  brief  in  the 
diffusive  regime.  This  would  not  apply  to  phase-trapped 
particles,  which  experience  an  extended  resonant  interaction 


Figure  4.  Bounce-averaged  quasi-linear  diffusion  coefficients  (solid  curves)  and  diffusion  coefficients 
for  coherent  interactions  with  monochromatic  waves,  averaged  over  the  same  frequency  and  wave  normal 
angle  distributions  (red  squares),  (top)  The  contributions  from  n  =  ±1  and  (bottom)  the  contributions  of 
just  n  =  0  are  shown.  As  predicted  analytically,  calculations  using  the  two  approaches  agree 
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time,  and  which  are  believed  to  be  key  for  the  self-consistent, 
nonlinear  growth  of  chorus  waves. 

[34]  For  computing  diffusion  coefficients,  there  is  no 
apparent  major  advantage  to  either  viewpoint,  the  same 
number  of  integrals  must  be  done  either  way.  However,  the 
coherent  interaction  approach  has  the  large  benefit  of  indi¬ 
cating  when  the  diffusion  approach  becomes  invalid,  and 
nonlinear  effects  must  be  considered.  Estimates  of  these 
effects  have  the  form  of  velocity  space  advection,  and  may 
be  included  in  a  combined  diffusion-advection  equation 
[Albert,  1993,  2000,  2002].  The  refinement  of  these  esti¬ 
mates,  and  their  use  in  global  simulations,  is  the  subject  of 
ongoing  work 


which  corresponds  to  equation  A8  of  Lyons,  and  which 
satisfies  equation  (Al)  above  for  any  choice  of  B 2(u;)  and 
In  the  notation  of  equation  (4), 


47tV 

utile 


Bi  avCGlG2> 


(A7) 


which  is  used  in  section  4. 
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Appendix  A:  Parameterization  of  the  Wave 
Distribution 

[35]  The  Fourier  transform  of  the  squared  wave  magnetic 
field  is 


02  =  f 

°wivc  / 


[Bk|2  dh. 
v  (2*)3 


(Al) 
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dl  k 
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{i*y 


(27r)3 
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1  duj 

0 

d0  kj , 
3 

for  any  cold  plasma  mode, 

d(k^k{) 


J  = 


d(u j.9) 


_ut  sm  6 


(A2) 


(A3) 


and  r(u;,  0)  is  /i2  /i  +  as  in  equation  (4) 

[36]  Lyons  [1974b]  explicitly  assumed  that  the  wave  dis¬ 
tribution  was  independent  of  both  <f>k  and  the  sign  of  so 
that  the  integrals  could  be  restricted  to  0  <  k\  <  oc  or  0  <  6  < 
7t/2,  with  an  additional  factor  of  2.  However,  to  connect  to 
single-wave  results,  it  is  more  natural  not  to  assume  sym¬ 
metry  with  respect  to  ±k  ,  and  to  take  0  integrals  from  0  to  tt. 
Then 


Following  Lyons ,  we  now  write 

bl  c=  r*  («)<**> 

Jo 

and  also  factor  B2(u;)  out  of  |Bk  2/V  This  leads  to 
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|Bk| 
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y  B2  (w  )  d  J 
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